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The Siberian forests are the habitat of many insect species. Periodic outbreaks of 
certain insect pests cause a decrease in tree growth, forest decline or mortality over vast 
areas. The Siberian silkmoth ( Dendrolimus superans sibiricus, Tschetw) feeds heavily 
on needles of certain tree species, rapidly defoliating and killing large stands of trees. 
Outbreaks seem to be influenced by topography. Remote sensing offers a tool for 
monitoring insect outbreaks over large areas however, not much has been reported in the 
literature about this. In this paper we looked at the use of data from the French SPOT- 
Vegetation satellite for insect pest monitoring. 

The spatial and temporal dynamics of an outbreak of the Siberian silkmoth were 
correlated with topographic features of the affected area using SPOT-VEG data and a 
high resolution digital elevation model (DEM). In 2002-2003 an outbreak affected 
-20,000 ha in the South Siberian mountains of Russia. The outbreak began between the 
elevations of -430- 480 m and on southwest slopes with steepness <5°. As the pest 
searched for food it moved up and down slope, resulting in an elevation distribution split 
within a range of -390-540 m and slope steepness up to 15°. In the final phase the 
azimuth distribution of damaged stands became even. The correlation between the initial 
phase and topographic features can be used to prioritize monitoring forest areas most 
vulnerable to destruction by pests. 

The significance of satellite-based pest monitoring should increase in the context 
of climate change. The nutrition base of Siberian silkmoth , so called “dark needle 
stands” composed by fir and Siberian pine, is within the zone of most intensive observed 
and predicted climate warming. Nowadays the zone of potential outbreaks is limited to 
-59 deg north latitude and by severe climate. Siberian silkmoth is a natural component of 



taiga ecosystems at higher latitudes too. As a response to warming, it’s expected that 
outbreaks will occur at higher latitudes, alongside with an increase of outbreak frequency 
and area of damaged stands in the south. Observed penetration of the Siberian silkmoth 
host species to the north, into currently larch dominated communities will, in the long 
term, lead to a vast increase a nutrition base of this pest. Due to the geographic scale 
and remoteness of taiga landscapes vulnerable to Siberian silkmoth outbreaks, Spot- 
VEG type sensors should be used in pest monitoring. 
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INTRODUCTION 

The Siberian forests are the habitat of many insect species. Periodic outbreaks of certain insect 
pests cause a decrease in growth increment, forest decline or mortality over vast areas. The Siberian 
silkmoth (Dendrolimus superans sibiricus, Tschetw) feeds heavily on needles of certain tree species, 
defoliating and killing large stands rapidly. This is one of the primary factors of taiga succession. 

The preferred pest host species are fir ( Abies sibirica ), Siberian pine ( Pinus sibirica ) but spruce 
(Picea obovata) and larch ( Larix sibirica) are also sometimes affected. Outbreaks are encouraged by 
favorable weather conditions: low summer precipitation, relatively mild winters with stable, dense 
snow cover and lack of late spring and early autumn frosts. On the contrary, cold, rainy summers and 
severe low-snow winters are not favorable for the Siberian silkmoth. Outbreaks have a periodicity, 
occurring about every 15-25 years. Between 1878 and 2004 ten Siberian silkmoth outbreaks were 
observed in the Yenisey River watershed area. The largest outbreak (1954-57) resulted in tree 
damage over about 4 million ha and tree mortality on about 1.5 million ha (Kharuk et al, 2003, 2004). 
The last outbreak occurred in the south Siberian mountain taiga in 2002-2003 yrs (fig. 1). Outbreaks of 
insect pests promotes wildfires, because pest-killed stands accumulate combustible material in form of 
dead wood, grass and bush communities. 
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Satellite data are an effective tool for detecting, mapping and monitoring of outbreaks (Nelson, 
1983; Dottavio and Williams, 1983; Royle and Lathrop, 1997, 2002; Gilli and Gorla, 1997). 

Spacebom instruments that provide higher resolution data (e.g., Landsat) are applicable for analysis of 
forest vigor and large-scale forest damage mapping (Luther et al, 1997; Radeloff et al., 1999; Kharuk 
et al., 2003). However, there are surprisingly few reports on coarse spatial resolution satellite data 
application for monitoring insect-induced tree defoliation and mortality (Gilli and Gorla, 1997; 

Kharuk et al., 2004; Fraser and Latifovic, 2005). The recent Siberian silkmoth outbreak (2002-2003) 
provides an opportunity to test the latest generation of coarse resolution spectroradiometers (SPOT- 
VEG, Terra/MODIS) for pest monitoring. Combination of this satellite data with available high- 
resolution 3D models (USGS, 2007) provides a tool for spatial analysis of outbreak dynamics which 
can reveal the relationship between outbreak patterns and relief features (elevation, azimuth, slope 
steepness). This can allow prioritization of monitoring for the detection of the initial phase of an 
insect outbreak,. Early detection is an important factor for effective pest control and outbreak 
suppression. 

The purpose of this paper is to analyze the spatial and temporal dynamics of a Siberian 
silkmoth outbreak based on SPOT-VEG data and a digital elevation model (DEM). 

Materials and method 

The most recent known Siberian silkmoth outbreak began in 2002 in the Eastern Sayan 
Mountains of Russia., resulting in tree mortality on 16,520 ha, according to on-ground data 
(Rosleszashita, 2005) (fig. 1). This area has severe continental climate. The mean annual temperature 
is -1.5 to -6,5 °C. Temperatures in the January, the coldest month, fall to -50°C while temperatures in 
July, the warmest month, may reach +17 to +19 °C. Mean annual precipitation is about 450-550 mm 
with the majority (60-70%) falling during summer. Mean relative humidity is 70-75% with maximum 
values in November-January and minimum in April-June. Stable snow cover lasts over 5 months 
normally from 1 November till 1 May. The growing season is 115-150 days. The outbreak was 
located in the vast watershed territory with relatively flat landscape with elevations up to 700 m. The 
dominant host species in the outbreak area were fir and Siberian pine. 
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Fig. 1. Siberian silkmoth outbreak location. Insert is SPOT-VEG based classification of 
outbreak area. 1. non-damaged stands, 2. outbreak area. 

Satellite (SPOT-VEG, Landsat) data, a forest inventory map (1: 100,000 dated by 1996 yr), and an 
insect damage map (of 1:100,000 dated by 2004 yr) were used in the analysis. The insect damage map 
illustrated two categories: insect-killed and non- damaged stands. Data were processed using Erdas 
Imagine and StatSoft Statistica software (Erdas Field Guide, 1999; StatSoft, Inc. 2003). The 
classification accuracy was evaluated by examing the classification confusion matrix and K-statistics 
(Erdas Field Guide, 1999). 

Image analysis 

SPOT-VEG data consisted of ten-day composites (24 scenes total, product S10, HDF format, 
Platte Carree projection, 16 bits radiometric resolution). The SPOT-VEG data were referenced to the 
middle of the ten-day compositing interval. Composites contained four spectral bands (430-470, 610- 
680, 780-790, and 1580-1750 nm) with 1 x 1 km pixel size. 

Ladsat-7 data (6 scenes) were used as an additional information source; this sensor is useful for 
host species and outbreak area detection (Kharuk et al., 2003). The analyzed scene fragments on the 
study area were essentially cloud free (< 5% of cloud). 

A high resolution DEM based on the US SRTM data was used for outbreak pattern analysis ( 
spatial resolution of 90 m and vertical resolution of ±15 m; USGS 2007). The outbreak area was 
analyzed with respect to eight slope aspect categories [(1) North (aspect angle range of 338°- 360° and 
0°— 22°); (2) Northeast (23°-68°), (3) East (69° - 1 13°); (4) Southeast (1 14° - 158°); (5) South (159° - 
- 203°); (6) Southwest (204°- 248°); (7) West (249° - 293°); (8) Northwest (294° - 337°)]. 

The distribution of topographic elements (altitude, azimuth and slope steepness) within the 
analyzed area is uneven; this could bias the results of the analysis. To make the analysis non-biased, 
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areas with a given damage category and with given azimuth, slope steepness and elevation were 
matched to areas with the similar parameters within the entire damaged area: 

^c(i) = A c(i)-f/ A. c(i)-I 

where the cq) subscript represents the ( h category of landcover feature c, A C (i)-f is the area of 
the given on-ground class within the i h category of the topography feature c, and A C (i)-i is the area of 
the t h category of topography feature c over the entire analyzed area. 

Data analysis included: 1) NDVI dynamics; 2) Damaged stand area dynamics; 3) correlation 
analysis of outbreak area dynamics and relief features (elevation, azimuth, slope steepness). Training 
sample selection was based on inventory and insect damage maps and Landsat data. The region 
growing method was used to form the training sample area. Since training sample statistics were 
found to be normally distributed, the maximum likelihood classification and threshold utility method 
were applied. Threshold utility allows to select study classes based on classifier distance and adjusting 
confidence level a (a = 0.01-0.3). Minimal calculated classifier distances for each pixel were used in 
analysis. 

The analyzed area included only two Siberian silkmoth host species (fir and Siberian pine). 
This allowed exclusion of background classes (such as deciduous stands, cuts, shrubs, etc), thus 
decreased classification error. For this purpose a “host species” mask was generated based on pre- 
outbreak SPOT-YEG scenes. 

The temporal set of SPOT-VEG based classification maps were generated that covered the 
period from 1999-2004 . Map generation included these steps: 1. Noise reduction by median filtering 
(3x3 window). 2. “Host species” stand mask application. 3. Removal from analysis previously 
damaged stands identified from existing insect damage maps. 

Table 1 and 2 represent the training area data and classification accuracy statistics. Note that for the 
beginning of the outbreak (Table 2, lines 2-4) the accuracy is poor . This is dues to the difficulty of 
classifying insect damage at the earliest stages. As the outbreak develops classification accuracy 
becomes satisfactory. (Table 2). 

Table 1. Ttraining set data for classification map generation 
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Table 2. Classification accuracy 
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5.13% 
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Damaged stands 
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0.92 
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Damaged stands 
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5.13% 
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Nf - number of reference points; N s - number of classified points; N c - number of correctly identified 
reference points; A c , A u - classifier and user accuracy; E 0 , E e - omission and commission errors; K c - 
current class K-statistics value; Ko - overall K-statistics value. 

Results 

The temporal dynamics of the NDVI and AND VI for pre-outbreak (1999-2001), outbreak 
(2002), and post-outbreak (2003-2004) periods is shown on fig. 2. 
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Fig. 2. NDVI dynamics of Siberian silkmoth outbreak and control areas. 

1 - ANDVI (a difference between control and outbreak stands), 

2 -NDVI of damaged stands, 3 - NDVI of control stands. 

Bars shows a confidence level (p<0.05). 

The NDVI profiles were generated based on temporal data from outbreak and areas of no 
insect damage (control areas) (Table 1). ANDVI is the difference between NDVI values for the 
control and outbreak areas. In three years prior to the outbreak (1999 - 2001) there were no significant 
NDVI difference between the outbreak and control areas. Therefore, on the time axis in Table 1 each 
year is represented by a single mid- August data point The ANDVI difference becomes significant 
starting end of July 2002 (fig. 2) when the insect outbreak began. In both the outbreak and control 
areas NDVI decreases in the mid August (08/15/2002) with a larger NDVI decrease in the outbreak 
area. This decrease is attributed to decrease in SPOT-VEG composite quality caused by haze on the 
ground. Starting mid-September the ANDVI differences stabilized. The NDVI decline observed at the 
end of September was caused by phenology changes. Note that ANDVI for the winter-acquired image 
(02/15/2003) is also high. During the following years’ growing seasons (2003 - 2004) the ANDVI 
decreased. 
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Fig. 3. Insect outbreak area over time 

The outbreak area plotted over time forms a logistic type curve with a lag-phase following by a 
rapid increase of defoliated area during next two years after which it levels off (Fig. 3). 

The elevation of the initial outbreak phase ranged from ~430 to 480 m as the outbreak 
continued, the elevation range expanded from ~390 to 540 m. (Fig. 4) During this period the area of 
the outbreak increment (the difference in outbreak area in relative units) the distribution becomes bi- 
modal. 
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Fig. 4. Change in outbreak area increment with respect to elevation. 

1, 2, 3 - increment of outbreak areas for the period of 08/1 1/02-08/21/02, 

08/21/02-09/01/02 and 09/01/02-09/11/02, correspondingly. 

The area increment was calculated as outbreak area distribution for given the dates. 

The insect outbreak began on gentle slopes of < 5°. As the outbreak expanded, areas with slope 
steepness up to 15° became affected (Fig. 5). The slope aspect of the outbreak area also changed over 
time. Initial damage was found on slopes with a southwestern aspect. Later, damage occurred more 
frequently on slopes with a western aspect. By the end of the outbreak, the slope aspect became more 
evenly distributed. (Fig. 6). 



Fig. 5. Outbreak area increment dynamics with respect to slope steepness. 

1, 2, 3 - increment of outbreak areas for the period of 08/1 1/02-08/21/02, 
08/21/02-09/01/02 and 09/01/02-09/11/02, respectively. 
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Fig. 6. Outbreak area dynamics with respect to slope aspect. 

Outbreak area 1 - 08/1 1/02, 2 - 09/01/02, 3 - 21/08/03 

Discussion 

In the area of study, a population increase of the Siberian silkmoth caterpillar was observed in 
1998 (Rosleszashita, 2005). However, it was not until 2002 that an outbreak of the insect pest 
occurred. Using SPOT-VEG data the first indication of Siberian silkmoth outbreak was detected at the 
end of July 2002, when the difference between NDVI of healthy and damaged stands became 
significant (Fig. 2). During the next three years of the active outbreak about -17 000 ha of forests 
were completely defoliated. A further increase in the outbreak area was due to mortality of weakened 
trees and, presumably secondary pest damage such as caused by the bark beetle ( Monochamus 
urussovi, Fish). In the years following the outbreak (2003, 2004) NDVI differences between outbreak 
and control areas decreased due to growth of grass, shrubs, and deciduous trees. . 

Because snow cover reduces the influence of the background the NDVI difference between 
healthy and damaged stands is also high during the winter. This indicates that winter acquired images 
are potentially useful in pest monitoring. NB: check and make sure this references a point in the prior 
text 

In previous insect outbreaks affecting the Sayan mountain region, the maximum insect induced 
stand defoliation is typically observed during late June through the beginning of July. The defoliation 
in this outbreak peaked later, during September vis a vis Figure 3 NBicheck fact. This time-shift could 
be attributed to two factors. First, the defoliation may have been caused by presence of more than one 
silkmoth generations (with one- and two-years life cycles) rather than one generation. Second, this 
outbreak occurred at relatively high elevations (390 - 520 m) compared to the majority of reported 
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outbreaks in the area which have been located in lowlands. Thus the impact of the local climate could 
be a factor in the timing of the peak of maximum defoliation. 

Earlier researchers empirically observed that the outbreaks begin in the best pest habitat sites 
(Isaev and Kondakov, 1986). The satellite data obtained shows that the spatial development of the 
outbreak was also related to topographic features. Because the initial outbreak phase occurred on 
gentle (<5°) southwest slopes within the elevation range of 420-500 m, these conditions appear to be 
the most favorable pest habitat (Fig. 4, 5). As food resources were depleted the pest migrated to less 
suitable areas with slopes up to 15° (Fig.5) and to both higher and lower elevations, expanding the 
elevation range to 390-540 m. This created a bi-modal distribution of elevation shown in Figure 4. 

In order to minimize forest damage, it is crucial to detect the initial stages of Siberian silkmoth 
outbreaks. The combined satellite-derived NDVI and DEM analysis revealed the topographic 
characteristics of sites where the outbreak began. This knowledge is essential for prioritizing areas to 
monitor for the first signs of silkmoth outbreak damage. It’s necessary to emphasize that the spatial 
pattern data are not biased, i.e. the normalization used provides results which are independent of the 
absolute areas with given topographic features. Moreover, data were normalized with respect to the 
area covered with host species. Otherwise the final result has to be dependant on the initial host species 
distribution: if this distribution is non- even with respect to topography, the resulting spatial outbreak 
pattern will be biased. 

The final SPOT-VEG derived outbreak area was about 20,000 ha, which is about 12% higher 
than the area reported by data collected in the field (Rosleszashita, 2005). This difference actually is 
within field collected data measurement error, which is aboutl0-15%. The outbreak area in this study 
is relatively small. The outbreak size was limited by weather conditions unfavorable for the insect pest 
species and, partially by suppression measures. In comparison, during a previous outbreak (1993- 
1995) that occured between the Yenisey and Angara rivers, the Siberian silkmoth affected stands on 
-700,000 ha and caused tree mortality on about 300,000 ha. This scale of damage allows applying 
NOAA/AVHRR data for outbreak analysis (Kharuk et al, 2004). The data analyzed here shows that 
even small-size outbreaks could be analyzed based on SPOT-VEG images: the outbreak area of about 
4 pixels in size (400 ha) was detectable (fig. 2, 3). Nevertheless, application of higher resolution data 
(e.g., Landsat, SPOT) is necessary for early detection of initial outbreak phases, which is important for 
effective pest control. Unfortunately high resolution data are non available with the needed frequency. 
There also are temporal “gaps” between data acquisition for the same area. For Landsat it is 16 days, 
the period comparable with active phase of outbreak dynamics (fig. 3). On the contrary, SPOT-VEG 
and Terra/MODIS data are available at the daily scale. 

An important problem is the indication of the intermediate levels of stand defoliation. In our 
study we were limited by maps derived from on-ground data with only two classes: healthy and 
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completely defoliated stands. Thus, only the dynamics of silkmoth-killed stands were analyzed. 
Whereas SPOT-VEG sensor should be sensitive to the less defoliated stand categories; this conclusion 
is based on the study of AVHRR application: also this sensor is less sensitive thans SPOT-VEG, three 
stand categories were detected (healthy, completely and intermediately defoliated; Kharuk et al, 2004). 
For the further progress it’ s necessary to combine timely organized on-ground studies with satellite- 
based observations during next Siberian silkmoth outbreak. 

Finally, the significance of satellite-based pest monitoring should increase in the context of 
climate change. The nutrition base of Siberian silkmoth , so called “dark needle stands” composed by 
fir and Siberian pine, is within the zone of most intensive observed and predicted climate warming 
(Zwiers, 2002). Nowadays the zone of potential outbreaks is limited to ~59 deg north latitude and by 
severe climate. Siberian silkmoth is a natural component of taiga ecosystems at higher latitudes too. 
As a response to wanning, it’s expected that outbreaks will occur at higher latitudes, alongside with an 
increase of outbreak frequency and area of damaged stands in the south. Observed penetration of the 
Siberian silkmoth host species to the north, into cunently larch dominated communities (Kharuk et 
al., 2005) will, in the long term, lead to a vast increase a nutrition base of this pest. Due to the 
geographic scale and remoteness of taiga landscapes vulnerable to Siberian silkmoth outbreaks, 
SpotVeg type sensors should be used in pest monitoring. 
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